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There is evidence for existence of massless Dirac quasi-particles in graphene, which satisfy Dirac 
equation in (1+2) dimensions near the so called Dirac points which lie at the corners at the 
graphene's brilluoin zone. It is shown that parity operator in (1+2) dimensions play an interesting 
role and can be used for defining conserved chiral currents [there is no 7* in (1+2) dimensions]. 
It is shown that the "anomalous" current induced by an external gauge field can be related to the 
anomalous divergence of an axial vector current which arises due to quantum radiative corrections 
provided by triangular loop Feynman diagrams in analogy with the corresponding axial anomaly in 
(1+3) dimensions. 

Recent progress in the experimental realization of a single layer problem of graphene has lead to extensive ex- 
ploration of electronic properties in this system. Experimental and theoretical studies have shown that the nature 
of quasiparticles in these two-dimensional system are very different from those of the conventional two-dimensional 
electron gas (2DEG) system realized in the semiconductor hctcrostructures. Graphene has a honeycomb lattice of 
carbon atoms. The quasiparticles in graphene have a band structure in which electron and hole bands touch at two 
points in the Brillouin zone. At these Dirac points, the quasiparticles obey the massless Dirac equation in (1+2) 
dimensions [l[. In other words, they behave as massless Dirac fermions leading to a linear dispersion relation = vk 
(with the characteristic velocity v ~ lO^m/s). This difference in the nature of the quasiparticles in graphene from 
conventional 2DEG has given rise to a host of new and unusual phenomena such as anomalous quantum Hall effects 
and a tt Berry phase P, Q • These transport experiments have shown results in agreement with the presence of Dirac 
fermions. The 2D Dirac-like spectrum was confirmed recently by cyclotron resonance measurements and also by angle 
resolved photoelectron spectroscopy (ARPEC) measurements in monolayer graphene Q. Recent theoretical work 
on graphene multilayer has also shown the existence of Dirac electrons with a linear energy spectrum in monolayer 
graphene [1]. 

The Dirac points lie at the corners of the graphene's Brillouin zone and have their position vectors in momentum 
space as @. 

K = |^(1,V3), K' = ^(l,-V3) 
da 6a 

where a is the carbon-carbon distance. 

Near the Dirac point K , the Dirac equation takes the covariant form Q 

z(7^9^)V^ - (1) 

where 

vf dt 

and 

7'^a^ = 7% + 7'ai+7252 (2) 

Now it is known Q that in 3 space-time dimensions there exists two inequivalent representations for 7-matrices 
[this is true for any odd number of space-time dimensions]: 

n SI ■ 1 9 ■ 9 

7^ = cr'^, 7 = ,T = 

7° = <j\ 7I =z(7\7' = -*<t2 (3) 
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One can take the second representation for the Dirac equation near the Dirac point K, which is obtained from K 
by the parity operation: 

Taking the two representations mentioned above into account one can write the parity conserving Lagrangian as 

L — i) j^(id)i)j^ + ■ip_{id)^- 

where 

d = 7"ao + 7'c>i +7^92 

d = -f"do + 7'9i - 7^92 (4) 
Parity operation takes the solutions in one representation to the other: 

r^{xP) = -Vp^+{x) (5) 
where = (a;°, , — x^). It is convenient to transform to new fields Q 

ipB = ij^ip- (6) 

The Lagrangian can then be written as 

C = i^Aih''d^)i;A+^Bih''d^)i^B (7) 
It is instructive to put mass term in the Lagrangian ([7]), which one can always put equal to zero: 

C = ij^{i-/''d^)^A + ■0b(*7^5^)'0b - mvfiip^ipA - ^b^b) (8) 
where under the parity operation 

Va,b{^^)^^W^bA^^) (9) 

The Hamiltonian density is 

= '^/[V'A(-«7'f^0^A + i'B{-'il^di)ipB + mvfiipA'^A - V'bV's)] (10) 

The Hamiltonian density has the so called conjugate symmetry Q, tpA ^ cr^i^B, in the sense that "H —H. 

It may be noted that the Lagrangian (j8]) is invariant, even in the presence of the mass term, two independent 
transformations 

i^A ^ e'"^VA, ^B ^ e^"-^B (11) 
where a a and ub are real, and has thus [/^(l) (g) [/^(l) symmetry. The corresponding conserved currents are 

J A = ^A^^A. J^B = i'Bl^'i^B (12) 

One can form even (odd) combination corresponding to "vector" ("axial vector") under parity 

J£ = 1/2[^^7^^A ± ^B^i^B] (13) 

In (1+2) dimensions, there is no 7^ available as in (1+3) dimensions. But still one may generate "chiral" currents. 
In fact under the infinitesimal transformations 

i^A,B -> -0^,3 = V'A.B + iaipB,A 

and 

^A,B iP'ar = ^A,B ± Olll)B,A 



the Lagrangian ([8]) respectively transforms into 



£2 = C-2mvf(^AipB + i'B'^A) (14) 



The corresponding conserved " chiral" currents in the absence of mass are 

1 



4 = T^i^An^B + ^Bri^A) (15) 
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1 

-2^ 



J5 = -^nii'Ari'B-i'BrM (16) 



These currents respectively correspond to "vector" and "axial vector" under parity. In the presence of mass term in 
the Lagrangian ([8|) 

df,j!^ = imvfitpj^tpB - i'Bi^A) (17) 

d^j!^ = mvf{ipj^i)B + i^B^A) (18) 

An external gauge field (electromagnetic) can be introduced by replacing the ordinary derivative by the covariant 
derivative 

-^D^ = d^ + -A^ (19) 

This gives the Dirac equation in (1+2) dimensions 

[i-f^'Df,Tmvf]^A.B = 0, (20) 
By multiply on the left by (^ij'^ D^, ^ rnvf) one can put the resulting equation in the Pauli form 

[D^'D^ + Yc'^'^'F^u + m'v}]ilJA, B - (21) 
This equation differs from the Klein Gorden equation in the term ^cr^'^ F^y, F/^^^ ~ d^A^ — d^A^. Using 

^ ^Am^^^ (22) 

one can write 



2c 



<J'"'F^, = 1/2/\a (23) 



where is current induced by the external gauge field A^ 



/A ^ ^e^^^'^F^, (24) 
c 

and has abnormal parity. The corresponding induced charge is 

Q = J d^x f{x) (25) 

where 

/O(^) = ^e^F = 2-(V X Ah - 2-B (26) 
c c c 

Here B is the magnetic field perpendicular to the x — y plane. Thus Q = where $ is the magnetic flux. 
The parity conserving Lagrangian which gives Eq. (21) is 

C± = [J^aID^'D,, + m\})^A - J^BiD^'D,, + m\})^jB] + l/2[(^ aI^'^a - i^Bl''M]U 

so that ffj, is coupled with the current given in Eq. (13) Q. The Lagrangian is invariant under J7a(1) ^ J7s(l). 
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Next we discuss whether can be related to "anomalous" divergence of some axial current, which arises due to 
quantum corrections. Obviously such a current can not be J^. Then the remaining axial current is given in Eq. 
(16b). Here analogy with the axial vector "anomalous" divergence in (1+3) dimensions, namely [lo| 



(9 



32n' 



pa 



(27) 



is useful. As is well known this divergence arises from quantum corrections provided by triangle graph which has two 
vector vertices and one axial vector vertex or that provided by its divergence 



where u and d denote up and down quarks (to is the quark mass), which provide the internal legs of the triangle. 
Note that although quark mass to appears above, but Eq. (27) is independent of quark mass. In our case we have 
corresponding ipA and ipB fields, which appear in the Lagrangian (8) or Hamiltonian (10). Noting that dfj,J^ involves 
mip a'Pb and mipgipAi the relevant Feynman graphs are shown in Fig 1 

fWWV 





FIG. 1: Triangle diagrams for "anomalous" current divergence. 

Note that it is essential to put mass transitions shown (see analogy with Majorana neutrinos) , such mass transitions 
are provided by the divergence of 9^ given in Eq. (17). Noting that A and B propagators involve opposite masses, 
the matrix elements are given by 

^ [ d^l , , I 1,1 1 , 

J (27r)-" /+ki-mvf /-mvf /- - tow/ /-k2+mvf 

+T~, — \ >(i — ? — ; — 7 )\\ 

/ + Ki + mvf / + mvf /— /C2 + mvf /— /:2 — mvf 

^ 2 2 2 [ Tr[{/- fix +mvf)j''{/ + mvf)'y'' ~ (/+ - mvf)■y^'{/ - mvf )Y] ki ^ k2 

f"" J {2^)D [{I + kif - m-^v}][P ^ m-^v}][{l - k2Y - m-^v}] ^ ^ u ^ 

The numerator in the integral which contributes is, 

iV^ = 277ivf Tr[YYl"{l + ki)p + YYl'lp] 
= -4mvfi[eP^''il + ki)p + e^P"lp] 

-Amvfie'^'Pkip (29) 

Using the Feynman parametrization, the denominator takes the the form, 

[{I + kix ~ k2yf - Af 

where 

A ~ m'^v'j — xk\ — yfcj + {xki — j/fc2)^ 

Making the shift I ^ I — [kix — k2y), the denominator becomes {P — A)^ and the dimensional regularization gives 
[D=3] 



T 



-4mVe^*.-'^''fc, f dx f dv ' r(3-i^/2) 1 , 



fci o k2 ^3Q^ 
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For photons on the mass shell 

kl = 0, fc2 = 

and puttmg 2ki.k2 — q^, A — m'^vj — q^xy. Thus we obtain, neglecting terms of order q^ 



Thus finally 



— (signm)e^''''(A:i-A;2)p (3f) 

ibTT 



^^,J^ = T7^(sign m)e^{k^)e,{k2)e^^''P{kl - k2)p (32) 

iOTT 



It may be noted while in (1+3) dimensions d^J^ is independent of to; in (1+2) dimensions, the corresponding 
quantity d^J^ is also independent of the magnitude of to but does depend on its sign (which is typical for odd 
space- time dimensions). In both cases mass was used as a regulator. In configuration space Eq. (32) takes the form 

= l^(sign ^)Axf^ = ^^(sign m)eP->^F^,Ay^ (33) 

[comparing with Eq. (27) in (1+3) dimensions]. Now in terms of electric and magnetic fields. 

e^'^'^F^.^Ax = -2A°B'^ - 2(1 x ~^)'^ = -2B[An + i^.^] 

where B is along z— axis and we have used ^ = x ^ so that = x yt. One can select a gauge in which 
^0=0 then Eq. (31) becomes 

= (34) 

[compare with second line (1^.1^) of Eq. (27) for axial anomaly in (1+3)]. It is instructive to also calculate dfj,J^ 
from the A— graph, the only change one has to make is to change the overall sign of the second term in Eq. (pS)) and 
m multiplying the integral to im. This changes N^- in Eq. p9|) to 



= + /ci)^r + (/ + kiYlf" -{1 + ki)l gP" + m^v]g^'''] (35) 
After making the shift I ^ I — {kix — k2y) and noting that terms linear in I do not contribute, one obtains 

T,. = I dxfj^ dy 5^''i^r(2 - d/2)(l)2-'^/2 + kt^k'( xy r(3 - d/2){^f-'^l^ 

k,^k2 ^ ^ 

where we have put kl = Q = fci.ei = fc2.e2=0 and then neglecting 2ki.k2 = q^, so that A = m'^v'j. Then finally 

11 

^M"^3 = (sign m){mvf 61.62 + (fc2.6i fci.62)} (37) 

47r 4 mv f 

Here the divergence does depend on to in addition to the (sign to). In terms of electric and magnetic fields given 
above (time independent A and uniform magnetic field) and neglecting the second term in Eq. ([37]). 

g2 

df,J^ = —(sign m)mvf B^r^ (38) 



One may take m = A/wj where A is energy gap in graphene's band structure, referred to as the Dirac gap 11 1 
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